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Abstract. We consider a system of PDEs of Monge-Kantorovich type arising from 
models in granular matter theory and in electrodynamics of hard superconductors. The 
existence of a solution of such system (in a regular open domain fl C R n ), whose con- 
struction is based on an asymmetric Minkowski distance from the boundary of Q, was 
already established in [G. Crasta and A. Malusa, The distance function from the bound- 
ary in a Minkowski space, to appear in Trans. Amer. Math. Soc.]. In this paper we prove 
that this solution is essentially unique. A fundamental tool in our analysis is a new 
regularity result for an elliptic nonlinear equation in divergence form, which is of some 
interest by itself. 



1. Introduction 

Let p: R n — ► R be a C 2 gauge function, i.e. a convex and positively 1-homogenous 
function, of class C 2 in W 1 \ {0}. In this paper we are concerned with the system of partial 
differential equations 

' -div(vDp(Du)) = f in ft, 
(1) < p{Du) < 1 in ft, 

p(Du) = 1 in {v > 0}, 

complemented with the conditions 

(u > 0, v > in ft, 
|ii = on <9ft. 

Here ft C W 1 is a bounded domain of class C 2 and / > is a bounded continuous function 
in ft. A solution of this system is a pair (u, v) of nonnegative functions, with u Lipschitz 
continuous in ft and v bounded and continuous in ft, satisfying the following additional 
conditions: (a) u = on 9ft; (b) p(Du) < 1 almost everywhere in ft; (c) u is a viscosity 
solution of p{Du) = 1 in the open set {v > 0}; (d) v is a solution of the first equation in 
(0) in the sense of distributions (see Definition |4.1| below) . 

This system of PDEs arises in some different situations. For example, the functions 
u and v can be interpreted respectively as the magnetic field and the power dissipation 
in a cylindrical hard superconductor of cross-section ft exposed to an external magnetic 
field linearly increasing in time (see e.g. ||). Moreover, in the case = |£|, (H)-© 
gives the stationary solutions of models in granular matter theory (see j|, ||). Another 
application concerns the existence of solutions to nonconvex minimum problems in calculus 
of variations (see J|, |?], |8| ) . Finally, Bouchitte and Buttazzo have studied a more general 
system in order to describe optimal solutions of some shape optimization problems. 



2000 Mathematics Subject Classification. 35C15. 

Key words and phrases. Distance function, Minkowski spaces, Hamilton- Jacobi equations, mass 
transport. 



2 



G. CRASTA AND A. MALUSA 



The results presented in this paper are an extension of those proved in 0, where the 
case p(£) = |£| was considered. An explicit solution to was constructed in |l(J. In 

order to describe this solution, we need some additional notation. Let cfo : CI — ► R be the 
Minkowski distance from the boundary dfl, defined by 

dn(x) = inf p°(x — y), x £ U, 

where p° is the polar function of p. It is well-known that dn is Lipschitz continuous in 
Q, and it is the unique viscosity solution of p{Dd^{) = 1 in Q vanishing on dfl (see fll6|| ). 
In it was shown that there exists a bounded continuous function vj: Cl — ► [0, +oo), 
whose explicit expression depends on /, p and on the geometry of Cl (see Section ||), such 
that the pair (dn,vj) is a solution to (|l|)-(|2|). 

The aim of this paper is to show that this solution is essentially unique. More precisely, 
we shall show that if (u, v) is a solution to ([l])-(^), then v = vj and u = dn in Qf = {x £ 
Q; Vf(x) > 0}. The proof of this uniqueness result is based on several ingredients. Some 
of them are an adaptation to our setting of arguments developed in [11, [§, 17, g|. A key 



point of the uniqueness proof consists in showing that, if (d n,v) is a solution to (H)-©, 
then v vanishes on the singular set S of dn (see Proposition |6^ below). In this respect, 
we use here a blow-up argument introduced by Evans and Gangbo for the case p(£) = |£| 
(see [11 1, Section 7, and ||), which in turn relies on the regularity of the solutions to the 
classical Laplace equation An = in an open set A C R n . In our setting the classical 
Laplace equation is replaced by 

(3) - div(Dp(Du)) = in A . 

Since the function a(£) := Dp(£) is defined and positively 0-homogeneous in R n \ {0}, no 
standard regularity result can be applied. In Section || we prove that, if p £ C 2 (R n \ {0}) 
and u is a Lipschitz continuous solution of (|3|) satisfying p(Du) = 1 almost everywhere in 



A, then u is of class C ,a locally in A (see Theorem 5A). Thanks to this regularity result, 



the blow-up argument of Evans and Gangbo still works in our setting (see Proposition 

2. Notation and Preliminaries 

2.1. Basic notation. The standard scalar product of two vectors x,y £ R n is denoted 
by {x, y), and \x\ denotes the Euclidean norm of x £ R n . Given two vectors v,w £ R n , 
the symbol v (8) w will denote their tensor product, i.e. the linear application from R n to 
R n defined by (v (g> w)(x) = v (w, x). 

By S" 1-1 we denote the set of unit vectors of R n , and by Mk the set of k x k square 
matrices. We shall denote by (ei, . . . , e n ) the standard basis of R n . The closed segment 
joining x £ R n to y £ R n will be denoted by [x,y], while (x,y) will denote the same 
segment without the endpoints. 

As customary, B t (xq) and B r (xo) are respectively the open and the closed ball centered 
at xq and with radius r > 0. 

Given A C R n , we shall denote by Lip (A), C(A), C b {A) and C k (A), k £ N the set 
of functions u: A — » R that are respectively Lipschitz continuous, continuous, bounded 
and continuous, and fc-times continuously differentiable in A. (Here and thereafter N will 
denote the set of nonnegative integers.) Moreover, C°°{A) will denote the set of functions 
of class C k (A) for every k £ N, while C k,a (A) will be the set of functions of class C k {A) 
with Holder continuous A;-th partial derivatives with exponent a £ [0, 1]. 

A bounded open set A C R n (or, equivalently, its closure A or its boundary OA) is of 
class C k , k £ N, if for every point xq £ dA there exists a ball B = B r (xo) and a one-to-one 
mapping ^: B -> D such that if) £ C k (B), ip- 1 £ C k (D), ip(B n A) C {x £ R n ; x n > 0}, 
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tp(B n OA) Cjiel"; x n = 0}. If the maps ip and ip~ l are of class C°° or C k > a {k G N, 
a£ [0, 1]), then A is said to be of class C°° or C ,a respectively. 

2.2. Convex geometry. By ICq we denote the class of nonempty, compact, convex sub- 
sets of IR n with the origin as an interior point. We shall briefly refer to the elements of Kq 
as convex bodies. The polar body of a convex body K G /Cg is defined by 

K° = {pG R n ; (p, x) < 1 Vx G if} . 

We recall that, if K G then K° G and K 00 = (K°)° = K (see |J, Thm. 1.6.1]). 
Given K G ICq we define its gauge function as 

p K (0 = mf{t > 0; £ G tK} . 

It is easily seen that 

Pk°(0 = sup{(£, p) ; p G K} , 
i.e. the gauge function of the polar set K° coincides with the support function of the set K. 
Let < ci < c 2 be such that ~B r -i(0) C if C £ _i(0). Upon observing that Upk(0 G if 

c 2 1 

for every £ / 0, we get 

(4) d|e| < PA -(0 <c 2 |e|, V£gM". 

We say that K G /Cq is of class Cl if <9if is of class C 2 and all the principal curvatures 
are strictly positive functions on dK. In this case, we define the i-th principal radius of 
curvature at x G dK as the reciprocal of the i-th principal curvature of dK at x. We 
remark that if K is of class C?, then if is also of class Cl (see |l^, p. 111]). 

Throughout the paper we shall assume that 

(5) K G /Cq 1 is of class C\ . 

Since K will be kept fixed, from now on we shall use the notation p = px and p° = p\. 

We collect here some known properties of p and p° that will be frequently used in the 
sequel. 

Theorem 2.1. Let K satisfy (||). Then the following hold: 

(i) The functions p and p° are convex, positively 1 -homogeneous in W 1 , and of class C 2 
in W 1 \ {0}. As a consequence, 

P (tO=t P (0, Dp(tH) = Dp(0, D 2 p{tt) = - t D 2 p{t), 

P°(tO = tp°(0, Dp°(tO = Dp°(0, D 2 p\tO = \D 2 p\i) ) 
for every £ G W 1 \ {0} and t > 0. Moreover 

(6) (Dp(0,0=P(0, (Dp°(0^) = P°(Z), V£GM"\{0} 

(7) D 2 P (0^ = 0, L» 2 p°(C)e = 0, V£G1T\{0}. 

(ii) For every £, ?7 G M n , we have 

(8) p(£ + »7)<p(f) + p07), P (e + ^)<P°(e) + P (r?) ; 

and equality holds if and only if £ and r\ belong to the same ray, that is, £ = A 77 or r\ = A £ 
for some A > 0. 

(Hi) The eigenvalues of the second differential D 2 p at v G S n ~ 1 are (with corresponding 
eigenvector v ) and the principal radii of curvature of dK at the unique point p G dK° at 
which v is attained as an outward normal vector. Symmetrically, the eigenvalues of D 2 p° 
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at v G S n 1 are ( with corresponding eigenvector v ) and the principal radii of curvature 
of dK at the unique point p G dK at which v is attained as an outward normal vector. 

3. Distance from the boundary 
Throughout the paper, we shall assume that 

(9) O C M. n is a nonempty, bounded, open connected set of class C 2 . 
Let us define the function 

(10) du{x) = inf p°(x-y), x eTt, 

that measures the distance from the boundary dtt to a point i £ O in the Minkowski 
norm associated to the polar function p° of p. Since dCl is a compact subset of M n and p° 
is a continuous function, the infimum in the definition of d^ is achieved. We shall denote 
by n(x) the set of projections of x onto dQ, that is 

(11) U(x) = {ye dn; d a (x) = p°(x - y)}, x G H. 

By abuse of notation, when H(x) = {xq} then we shall use II(x) to indicate the point xq. 
It is well-known that dn is a viscosity solution of the Hamilton-Jacobi equation 

(12) p(Du) = 1 in O. 



More precisely, it is the unique viscosity solution of (12) satisfying the boundary condition 
u = on dVL. Moreover, dn(x) > for every x G 0, dn G Lip(J7), and p{Du{x)) = 1 for 
a.e. (see |, H). 

We say that x G SI is a regular point of £1 if H(x) is a singleton. We say that x G £1 
is a singular point of if x is not a regular point. We denote by S C S] the set of all 
singular points of 0. It is well-known that cfo is differentiable at every regular point of £1 
(see ||, ||, 16]; see also Theorem |3.1|(i) below). 



From now on, for every xq G d£l we shall denote by ki(xq), . . . , K n _i(xo) the principal 
curvatures of 90 at x$, and by v{xq) the inward normal unit vector to dQ at Xq. We 
extend these functions to O \ S by setting 

= ^(II(x)), Kj(x) = Kj(II(cc)), i = 1, . . . ,n — 1, xGil\S. 



We collect in the following theorem all the results proved in [10] that are relevant to 
the subsequent analysis. 

Theorem 3.1. Let and K satisfy respectively @ and (||). Then the following hold. 

(i) Ecfi, and the Lebesgue measure o/S is zero. 

(ii) Let x G O and xo G n(x). Then, for every z G [xq,x), du is differentiable at z and 
(13) Ddh(z ) = J^2) . 

p(n^o)) 

f m j T/ie function dfi is of class C 2 onfl\S. 

Proof. See Remark 4.16, Corollary 6.9, Lemma 4.3 and Theorem 6.10 in [|n]]. □ 

At any point xq G d£l there is a unique inward "normal" direction p(xq) with the 
properties IL(xo + tp(xo)) = {xo} and dfi(xo + tp(xo)) = t for t > small enough (see p^| , 
Remark 4.5]). More precisely, these properties hold true for t G [0, r(xo)), where t(xq) 
is the normal distance to the cut locus S, defined below (see |i~0| , Proposition 4.8]). It 
can be proved that p(xq) = Dp(u(xo)) (see |l^, Lemma 2.2] and [JToL Proposition 4.4]). 
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From Theorem |3.1| (ii) and the positive O-homogeneity of Dp it is clear that Dp(v(xo)) 
Dp(Ddci(xo)). Summarizing, given x G O we have that 

(14) xq G IL(x) <^=^ x = xo + t Dp(Ddfi(xo)) for some t G [0, t(xq)], 

and, in such case, cfo(x) = 

The above considerations motivate the following definition. 



Definition 3.2. The normal distance to cut locus of a point x G O is defined by 
(15) t(x) 



mm{t > 0; x + tDp(Dd n (x)) G £}, if x G fi \ E, 
0, if x G E. 



The cut point m{x) of x G f2 \ £ is defined by m(x) = x + T(x)Dp(Dd^i(x)). 

Proposition 3.3. Let Q satisfy @. T/ien r zs continuous in fL Furthermore, there exists 
p > suc/i i/iai t(xo) > /or every xo G <9$7. 

Proof. Sec pi, Lemma 4.1 and Theorem 6.7. □ 



From Theorem |3 . l| (iii) , the function dn is of class C 2 on £2 \ E. We can then define the 
function 

(16) W(x) = -D 2 p(Dd n (x))D 2 d n (x), xgH\E. 

For any xo G d£l, let denote the tangent space to d£l at xo- If x G O \ E and LI(x) = 
{xo}, we set T x = T XQ . Observe that, by (|l3|) and @, we have D 2 p{Dd^i{x))v{x) = 0. 
Then, for every v G T x , we have W{x) v G T x . Hence, we can define the map 

(17) W(x):T x ^T x , W(x)w = W(x)w, 

that can be identified with a linear application from R n_1 to M n_1 . 
We shall use the following results (see [[lG|], Lemmas 4.10 and 5.1). 

Lemma 3.4. Let xo G dil. Then 

det[J -tW(xo)] = det[/ n _i - tW(x )\ 

for every t G M., and both determinants are strictly positive for every t G [0, t(xq)). 

Remark 3.5. Let xo G d£l. We recall that Dp(u(xo)) is an eigenvector of W{xq) with 
corresponding eigenvalue zero (see [|l(], Lemma 4.18]). On the other hand, from Lemma |3~4] 
we deduce that a number k ^ is an eigenvalue of W(xq) if and only if it is an eigenvalue 
ofPF(xo). 

Although the matrix W{xq) is not in general symmetric, its eigenvalues are real numbers, 
and so its eigenvectors are real (see [10], Remark 5.3). The eigenvalues of W(xq) have an 
important geometric interpretation. 

Definition 3.6 (p-curvatures) . Let xo G dVt. The principal p-curvatures of dQ at xo, 
with respect to the Minkowski norm dn, are the eigenvalues ki(xo) < • • • < k n -i(xo) of 
W(xq). The corresponding eigenvectors are the principal p-directions of dfl at xq- 



Up to now we have analyzed some properties of the matrices VF(xo) and W(xq) at points 
xo G dO,. Now we are interested in the evolution of these matrices along the transport ray 
starting from xq- 
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Lemma 3.7. Let xq G dfl, and define 

(18) V(t) = W(x + tDp(Dd n (x ))), V(t) = W(x + tDp(Dd n (x ))) 
fort G [0, r(xo)). Then 

(19) V(t) [I - t V(0)] = V(0), V(t) [/„_! - tV(0)} = 7(0), 

/or even/ i G [0,r(xo)). Furthermore, TrV(t) = TrV(t) for every t G [0, r(xo)). 

Proof. Let us consider the principal coordinate system at xo, i.e. the coordinate system 
such that xq = 0, e n = z^(xo) and coincides with the i-th principal direction of <90 at 
Xq, i = 1, . . . ,n — 1. Let X : U ^ M. n be a local parametrization of 50 in a neighborhood 
of = 0. The relation (|i~3|) can be written an 



(20) Dd a {X{y) + 1 Dp(Dd n (X(y)))) = ^ , y£W, t G [0, r(X(y))), 

where iV(y) = Recall that 

<9iV 

(21) N(0) = e n , — (0 ) = -«iei, i = l,...,ra-l, 

<% 

where kj, . . . , K n _i are the principal curvatures of 90 at xq. Differentiating ( p0|) with 
respect to y, at y = and using ( ^Tl) we obtain, for every i = 1, . . . , n — 1, 

(22) D 2 dn(x + tDp(Ddn(xo))) [i + tD 2 p(Dd Q (x ))D 2 dn(x ))] e; = 

where £j = [p(^) — (Dp(u), a) v]/p(v) 2 , and ^ = i/(xo). Differentiating (20) with 
respect to t at y = we get 

(23) D 2 d n (x + t Dp(Dd u (x ))) Dp(Dd u (x )) = . 

Let us apply D 2 p{Dd$i{xQ + t D p(Ddn(xo)))) = D 2 p(Ddo,{xo)) to both sides of equations 
(|22l ) and (|2^). Recalling the definition (|l6| ) of PF, we obtain the relations 

W{x + tDp(Dd Q (x ))) [I-tW(x )] = i = l,...,n-l, 

W r (xo + t-Dp(Ddn(a:o)))£>p(£)do(a;o)) = 0, 

where ^ = D 2 p(Ddn(xo)) £j. We have that, for every i G [0, r(xo)), 

(V(t)[I-tV(0)]ej = i = l,...,n-l > 
\y(t) [/-tF(0)] Dp(Dd n (a:o)) = > 

where in the second identity we have used the fact that V(0) Dp(Dda(xo)) = (see |lC| , 
Lemma 4.18]). Since (Dp(Ddn(xo)), u) = p(u) > from the positive 1-homogeneity of p, 
it follows that the vectors e±, . . . , e n _x, Dp(Ddn{xo)) span M. n , hence from (p4[) we infer 
that the matrix V(t) [I — t V(0)] is independent of t. Evaluating this matrix at t = 0, we 
finally obtain the first identity in ((H|). From ([l^) and (0) we have that 

(V(t) v, v) = -p(u) (D 2 d n (x + tDp(Dd n (x ))) v, D 2 p(v) u) = 



(24) 



for every v G M n , hence the second identity in ( )19| ) is also satisfied, and Tr V(i) = TrV(t). 

□ 



Remark 3.8. Let V(i) be the function defined in (18). By definition, the eigenvalues 
, re n _i of V(0) are the principal p-curvatures of <90 at xq, and the corresponding 
eigenvectors w\, . . . ,w n -i are the principal p-directions of <90 at xq. From the identity 
(|l9l) we obtain that 

V(t) (1 - i Kj)^ = KiWi, i = 1, . . . , n - 1, t G [0, r(x )). 
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Since for every t £ [0, t(xq)) the point xo+t Dp(DoIq(xq)) belongs to f2\£, then 1— t iki > 



sec 



10, Lemma 5.4]), and hence the eigenvalues of V(t) are 

K i 



ki(t) 



1 — tiki ' 



i = 1, . . . , n — 1, 



with corresponding eigenvectors w\, . . . , w n -\. 

Proposition 3.9. For every xq G 50 Zei its define the function 



(25) M so (a,t) = exp(- / TrW(x + a Dp(Dd n (x )) da ) , 

/or s,i £ [0, r(a;o))- TVien 



(26) 



M , t) = detW(x + sDp(v(xo))) = det 
detW(x + tDp(u(x ))) det 



/ n _i - w(x ; 



n-1 



n-1 



■pr 1 -tikj(x ) 

M 1 - SKj(x ) 



i=l 



for every s,t G [0,r(xo)) ; where iki(xo),--- ,ik n -i(xo) are the principal p-curvatures of d£l 
at Xq. 



Proof. Let V(f), i G [0, r(xo)), be the function defined in fllq). Prom Lemmas 3^4 and 3/7 
we have 

(27) F(t) = F(0) [I„_i - tV(0)]-\ Vt G [0, r(x )) . 

This implies that the matrix- valued function V(i) satisfies the differential equation 

V'(t) = V(t) 2 , \/t G [0, t(xq)) , 

and hence the function = det V(t) is a solution of the differential equation S'(t) = 
Tr[V(t)] 5(t) in [0,t(xq)). Now, the first equality in ( |26"1) follows from the fact that 

6(t) = 8(s) exp ( [ TrV (a) da) = 5(s)/M X0 (s,t) , 



for every s,t £ [0,t(xq)). The second equality follows from (p7|) . The last equality is a 
direct consequence of the fact that iki(xo), . . . , ik n -i(xo) are the eigenvalues of the matrix 

F(o). □ 



Proposition 3.10. The function M Xo (s,t), defined in §23), is jointly continuous with 
respect to xq G <90 and s,t G [0,r(xo))- Furthermore, 



n-1 



(28) 

where 

(29) 

(30) 



< M X0 (s,t) < Y[(1 + TK-) , Vx G 90, < s < t < t(x ), 



i=l 



T = max{r(x); x G <9f2}, 
if_ = max{[Kj(a;)]_; x £ dft, i = 1, . . . , n — 1}, 



6em<7 [a]_ = max{0, —a} i/te negative part of a real number a. 
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Proof. The continuity of M Xo (s, t) follows from its definition and the continuity of Kj. The 
estimate ( p8[ ) follows from the representation formula (^6|) and the estimate 

that holds for every xo £ dfl and < s < t < t(xq). □ 
For every x £ fl \ S let us define the function 

(31) M x (t) = exp (- j TtW{x + a Dp(Dd Q (x)) da^j , t £ [0, t(x)) , 



where W is the matrix defined in (|17|). For an explicit computation of M x it can be of 
some aid to recall that 

TrW(x) = TrW{x) = - Tr \d 2 p(Dd n (x)) D 2 d n (x 

for every x £ f2 \ £ (see Lemma 3.7). Given x £ O \ S, let 11(2;) = {xo}. By ( |l4| ) we have 
that 

W(x + a Dp(Dd Q (x)) = W{x + (cfo(x) + a) Dp(Dd n (x)), a £ [0, r(x)), 
which implies the relation 

(32) M x {t) = M X0 {d n {x),d n {x) + t), Vt£[0,r(x)). 

From the identity (p6|) we have that 



(33 ) A4w = n 1 '' rf " ( 1|' ) A (I) 

where Kj(x) := Kj(xo), i = 1, . . . , n — 1, are the /^-curvatures of dfl at xo- 

4. Existence of solutions 

In this section we recall the existence result for system (0)-(§) proved in Thm. 7.2]. 
The rigorous meaning of solution is the following. 

Definition 4.1. Let ft C M. n satisfy @ and let / £ Cb(fl) be a nonnegative function. A 
solution of system (|])-(|]) is a pair (u,v) of functions satisfying the following properties: 

1. u £ Lip(J7), u £ Cft(r2), u, w > in Vl; 

2. u = on dfl, p(Du) < 1 a.e. in fl, and u is a viscosity solution of 

p(Du) = 1 in {u > 0} ; 

3. u is a solution in the sense of distributions of — div(v Dp(Du)) = f in ft, that is 



(34) / v(x) (Dp(Du(x)), D(p(x)) dx = / f(x)<p(x)dx 

Jo, Jn 

for every p £ C£°(Q). 

Remark 4.2. Since v £ L°°(Q) and Dp(Du) £ [L°°(J7)] n , by a standard density argument 
|) holds for every test function ip in the Sobolev space W ' (fl). 

For / £ Cfe(fi), let us define the function 

r r(x) 



(35) v f (x) 



f(x + tDp(Dd n (x)))M x (t)dt, if x £ fi\E, 

if x £ E . 
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Theorem 4.3 (Existence). Let Q C W 1 satisfy (@) and let f > be a bounded continuous 
function in Q. Then, the pair (dn, vt) is a solution to (|l])-(||) in the sense of Definition 4-1 



A rigorous proof of Theorem |4.3| was given in [jlO| , Theorem 7.2]. In Section y we shall 
prove that the pair (du,Vf) is essentially the unique solution to problem In order 



to gain some insight in the representation formula (|35|), a formal derivation of (35) might 
be in order. 

Assume that (dn,v) is a solution of (H)-©, and that v G C 1 (r2 \ E), with v vanishing 
on E. Outside E, the equation — div(v Dp(Du)) = f is satisfied pointwise, that is 

u(a;) TrW(x) - (Dv(x), Dp{Dd n {x))} = f{x) , 

where W(x) = —D 2 p{Dd$i{x)) D 2 dn(x). Furthermore, from Lemma 3.7 we have that 
TrW(x) = TrW(x), where W(x) is the matrix defined in (|T 



Let x G Q \ E, and define = v(x + 1 Dp(Ddn(x))), t G [0, r(x)]. The function f(i) 
satisfies the following linear differential equation 

v'(t) = [TrW(x + tDp(Ddn(x)))] v(t) - f(x + t Dp(Dd n (x))) 

in [0, t(x)], supplemented by the boundary condition 

v{t(x)) = 

since x + t{x) Dp{Ddo J {x)) G E. The solution of this Cauchy problem, evaluated at t = 0, 
gives 

v ( x ) = v(0)= [ TX f{x + tDp(Dd Q (x)))M x (t)dt, 



o 

that is, the solution v (x) has to be the function defined in formula (|35|). 

With this heuristic in mind, our aim will be to prove that, if (u,v) is a solution to 
(jl])-(||), then (dQ,v) is a solution too (see Lemma |6.3|(ii)), and that if (cfo,i>) is a solution 
to (|l|)-(|2|), then v must vanish on E (see Proposition |6.7| ). The first goal will be achieved 
using the same arguments proposed in ||, whereas the second one needs a new regularity 
result for solutions of elliptic equations, which seems to be of some interest by itself, and 
that will be proved in the following section. 

5. A REGULARITY RESULT 

The aim of this section is to prove the following regularity result. 

Theorem 5.1. Assume that p is the gauge function of a convex body K satisfying 
Let A C W 1 be an open bounded set, and let u G W 1,oc {A) be a solution in the sense of 
distributions of the equation 

(36) - div(Dp(Du)) = in A, 
that is 

{ (Dp(Du(x)), Dip(x)) dx = 
J A 

for every (p G C£°(A). If in addition 

(37) p(Du(x)) = 1 a.e. in A, 
then u G Cj£(A). 

We recall the standard regularity result about solutions to the equation 

(38) - div(a(L>u)) = in A, 
(see pi §4.6], pi, §8.2]). 
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Theorem 5.2. Assume that the vector-valued function a belongs to C 1 (A) and satisfies 
the following growth conditions: there exist p > 1, and ao,/?o > such that for every 

(39) KOI + (1 + \Z\ 2 )*\Da®\ < a (l + |£| 2 )^ , 

(40) (Da(£)w, w) > foil + leP^M 2 , VweE". 

T/ien every solution u G W 1 ' 00 ^) o/ @) 6don 5 s to Cj£(A). 

In our case, p is a positively 1-homogeneous function of class C 2 in W 1 \ {0} satisfying 
the bounds @, and then Dp is a positively 0-homogeneous function of class C 1 in K n \{0}, 
but in general p is not even differ entiable at the origin. Moreover, ([!(]) cannot be verified 
near the origin. Hence we have no chance to apply directly the standard regularity results 
to solutions of equation (|36|) . 

An easy trick in order to have the "right" growth is to consider the function 

(41) l{t) = \p{t?, (GK". 



We have that 

(42) D 7 (0 

Clearly, a function u G W 1,co (A) satisfying (|37]) is a solution to (^) if and only if it is 
a solution to 



0, if£ = 0. 



(43) 



- div(Dy(Du)) = in A, 
p(Du) = 1 a.e. in A. 



Moreover, 7 (£) is a positively 2-homogeneous function of class C 2 in lR n \ {0}, hence 
D 2 j(£) is 0-homogeneous and continuous in 1" \ {0}. In particular the matrix-valued 
function D 2 ^ is bounded in W l \ {0}, and 

(44) p 2 7 || < c 5 = . max max{ |Aj7(0li £ ^ S"" 1 } . 

The following positive constants will be used throughout this section: 

c 3 = max{\Dp(0\; (£R"\ {0}} = max{|Dp(i/)|; 1/ G S^ 1 }, 
c 4 = max{\D 2 p(u)\; ^S"- 1 }, 

(45) 

r = min{r^(p); p G , i = 1, . . . , n — 1}, 
R° = max{rf (p); p G dK°, i = 1, . . . , n - 1} 

where rj(p) < . . . < r^_ x (p) are the principal curvatures of dK° at p. We remark that 
r° > since if is of class C\. 

The first technical tool is to prove that 7 satisfies some growth conditions similar to 
PD and fl4|) with p = 2. 



Lemma 5.3. Lei 7 6e i/ie function defined in (f4l|). T/ien 

(46) |Z) 7 (0| <c 2 c 3 |e|, 

(47) <^7(£),0>c 2 |£| 2 , 
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for every £ 7^ 0. Here c\, c%, C3 are the constants defined in (Q) and (|4q). Moreover there 
exists a constant cq > 0, independent of £ 7^ 0, suc/i i/iai 

(48) (.D 2 7(0^, > c 6 \w\ 2 , 
for every w G M n . 

Proof. Since 7 is 2-homogeneous, by Euler's formula {Dj(£), £) = 2j(^), and then 

(D^),0=p(0 2 >cl\C\ 2 . 

On the other hand, by (|42|), 

|£> 7 (oi <kowoi <Q2C3iei- 

It remains to prove that there exists a constant cq > such that 

(49) ^L> 2 7(£)w, w\ > c 6 \w\ 2 , Vu>GlR n , V£^0. 

Since 7 is a convex function, of class C 2 in M n \ {0}, we have that the quadratic form 
w 1 — > (D 2 ^/(^)w, w) is positive semidefinite for every £ 7^ 0. We shall show that, in fact, it 
is positive definite uniformly with respect to £ 7^ 0. 
Fixed £ 7^ and using the notation f = £/|£|, we have 

L> 2 7(£) = £> 2 7(^) = piy) D 2 p{v) + Dp(i^) ® Dp(u), 

so that 

^Z? 2 7(C)u>, w) = (d 2 j(v)w, = p(u) ^D 2 p(v)w, w\ + (Dp(u), w) 2 
for every w G M n . 

Fixed w G R n , let us denote by A = (v, w), and by w the projection of w on the 
orthogonal space C^{v) to ^, so that w = w + Xu, (w, v) = 0, and \ w\ 2 = \w\ 2 + A 2 . 

From Theorem [2.1| we have that D 2 p{v)v = and D 2 p(v) is positive definite in C^~{v) 
with 

(d 2 P (v)w, uj) > r°|u7| 2 , 
for every 117 G £^(z^), where r° > is the constant defined in p5|). Hence we obtain 

(^D 2 p(v)w, w\ = ^D 2 p(u)w, > r°\w\ 2 , 
for every w G 1™. On the other hand, by the 1-homogeneity of p we have 

(Dp(u),w + \v) = (Dp(u), w) + \p{v) . 
Hence we get the inequality 

(D 2 j(0w, w) > r°p(v)\w\ 2 + (Dp(v), w} 2 + \ 2 p{v) 2 + 2\p(v) (Dp(v), w) , 

for every w = w + Xu G M n , w G C^{v). 

It remains to prove that there exists cq > independent of v such that 

r°p(u)\w\ 2 + (Dp(u), w) 2 + \ 2 p{v) 2 + 2\p{u) (Dp(u), w) > c 6 (\w\ 2 + A 2 ) . 

We have that 

ypHM 2 + (D P (U), W) 2 > (l + y r^rp) 

(50) " ; \ ' 



,2 



12 G. CRASTA AND A. MALUSA 

where c\ > and C3 > are defined respectively in (||) and fl45"l) . If c = 1 + j from 
we get 

r°p{v)\w\ 2 + (Dp(i/), HI) 2 + \ 2 p{vf + 2Ap(i/) <Dp(i/), w) 



> 
~ 2 



^(-)|W| 2 + c (W), TO) + ^) + (l - ^) A 2 p(-) 
>^M 2 +(l-i) c 2 A 2 > C6 (|w| 2 + A 2 ) , 



where c 6 = min (l - A c 2 ) > 0. □ 

Thanks to the estimates (^6|), (|47|), and (^) we can now prove that a solution u to the 
first equation in (ft3| ) belongs to ^ 2 C (A). This part of the proof is based on a standard 
argument in regularity theory (see e.g. [|l^, §8.2]). The only point that should be stressed 
concerns the regularity of 7. Namely, 7 is of class C 2 in W 1 \ {0}, but in general the 
positively 0-homogeneous function D 2 7 is not even defined at the origin. In our case 
this is not a real problem, since the condition p{Du{x)) = 1 guarantees that Du(x) stays 
always outside a ball centered at the origin. 

Lemma 5.4. Every solution u £ W 1,co (A) of fcffi ) belongs to H 2 oc (A). 
Proof. The function u solves 

(51) / (Dj(Du(x)), Dtp{x)) dx = 



for every cp £ C£°(A). Let us fix <p, a coordinate direction e^, and /i 6 1 such that 
ifh{x) = <p( x ~~ h e i) has compact support in A. Choosing tp^ as a test function in (51) and 
making a change of variables in the integral we obtain 



(52) / (Dj(Du(x + heA), Du>(x)) dx = . 

J A 

Taking the difference between (|5l| ) and (^), we get 

(53) / (Dj(Du(x + hei)) - Dy(Du(x)), D<p(x)) dx = . 



On the other hand we have that p(Du(x)) = p(Du(x + hei)) = 1 for a.e. x G A. Then 
the function a(t) = D^{{1 — t)Du{x) + tDu(x + hei)) is continuous in the interval [0, 1]. 
In addition, either a S C 1 ([0, 1]), or there exists to £ (0, 1) such that a is of class C 1 and 
with bounded derivatives in [0, 1] \ {to}. Hence a(t) is a Lipschitz function, and 

D~f(Du(x + he)) - D-y(Du(x)) 



(54) 



r 1 d 

/ -r[Dj({l-t)Du(x)+tDu(x + hei))}dt 
Jo dt 

[ D 2 j((l - t)Du(x) + tDu(x + hei))[Du(x + ha)) - Du(x))]dt, 
Jo 



for a.e. x S A. Finally, the matrix 



Lh(x) 



[ D 2 j((l-t)Du(x)+tDu(x + hei))dt 
Jo 



is well defined and, for w € W 1 , the integrals of the kind 

(L h (x)w, w) = j (-D 2 7((l - t)Du(x) + tDu(x + ha))w, in) dt 
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satisfy for a.e. x G A 

(55) c 6 \w\ 2 < (L h (x)w, w) , Vwe R n , 

(56) ||L h (x)||<CB, 

where cq and C5 are the positive constants defined in Lemma |5.3| in (|44|). 

If we denote by A/ l n(x) = HiE+feglzliEi ; by ( [54] ) equation ( |53[ ) can be rewritten as 



(57) / (L h (x)DA h u(x), Dip(x)) dx = 

JA 

for every ^ G C£°(.A) and for every \h\ < dist(suppy?, cL4). By a density argument, 
we have that (|57| ) remains valid for every test function in Hq(A'), where A' is an open 
set compactly contained in A, and for every \h\ < dist(A', dA). Hence we can choose 
ip = ?y 2 Afrit as test function in (|57|), where r\ G C£°(A) is defined in the following way: 
given xo G A and r > such that the ball B2 T {xq) is compactly contained in A, we require 
that < r)(x) < 1 in A, rj(x) = in A \ B2 r (xo), rj(x) = 1 in B r (xo), and there exists 
m > such that \Drj(x)\ < m/r in A. With this choice of the test function, (|57|) becomes 



(58) 



(L h (x)DA h u(x), DA h u(x)) r? (x) dx 

(L h (x)DA h u(x), Drj(x)) ■q{x)A h u{x) dx , 
for every < dist(xo, cL4) — 2r. Recalling ( |55| ) we have that 
(59) c 6 / |DA /l u(x)| 2 r/ 2 dx < / (L h (x)DA h u(x), DA h u(x)) rj 2 (x) dx 



On the other hand, by (|5q ) and Young's inequality we obtain that there exists a constant 
C > such that 



(60) 



(L h (x)DA h u(x), Dt](x)) r](x)A h u(x)dx 

< C (tJ \DA h u(x)\ 2 r] 2 (x)dx+ ^ | J Dr ? (a;)| 2 |Afrn(x)| 2 ^ , 

for every e > 0. Choosing e small enough, from (|58|), (|5G|), (|60|), and the estimates of |£>f?|, 
we get 

M 



(61) / \DA h u{x)\ 2 dx<-^ [ \A h u(x)\- <I.r . 

JB r (x ) H JB 2r (x ) 



which implies, by a standard argument, that u G Hf oc (A) (see Lemmas 7.23 and 7.24 in 
@, or || §8.1] ). ' □ 



Proof of Theorem 5A. We already know that a solution u G TV ,00 (^4) of (|371)-(|36f) is also 

d.r, 



a solution to (f43|). Fixed 93 G C^°(r2), we can choose Jjj^ as test function in the weak 



formulation (pl|), obtaining 

(62) / IdMDu(x)), D^r-(x)\ dx = 

J A \ OXi I 

Moreover, by Lemma [TJ, the function u belongs to Hf oc (A). Then an integration by parts 
leads 

(63) J ^D 2 -f(Du{x))D-^-(x), Dip(x)\ dx = . 
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for every <p G C£°(Q). Then the partial derivative is a bounded solution of the linear 
elliptic equation 

(64) - div(D 2 j(Du(x))Dv) = in A, 

where the matrix D 2 j(Du(x)) satisfies the hypothesis of the De Giorgi-Nash regularity 
result (see pi §3.14], Theorem 14.1, pfl, Theorem 8.22). Hence we can conclude that 



the partial derivatives of u are locally Holder continuous. □ 

Remark 5.5. Let us define the function 7 P (£) = p{Cj p /Pi V > 1- Clearly a function 
u G W 1 '°°{A) is a solution to (|37|)-(|36|) if and only if it is a solution to 

f-div(D7 P (D«(z)))=0 in A, 
1 J |/)(Du(x)) = 1 in A 

where, as usual, the first equation is interpreted in the sense of distributions and the 



second in viscosity sense. Arguing as in the proof of Lemma 5^3 , it can be proved that 
there exists a positive constant c p such that, for every £ ^ 0, 

\d 1p (o\< ^ler 1 , 

(d 1p (0,0>^\ p , 

'D\(0 W > w ) > Cp\£\ P ~ 2 M 2 , G R n . 



We remark that, if p > 2, then j p is of class C 2 on W 1 , and the estimates above hold 
for every £ G R n . On the other hand, if p > 2, we cannot obtain an estimate of the 
type ( f40| ) near the origin, due to the p-homogeneity of the function j p , which implies the 
(p — 2) -homogeneity of -D 2 7 p (£). For this reason we have considered the case p = 2. 

6. Uniqueness 

This section is devoted to the proof of the following uniqueness result. 
Theorem 6.1. Let (u,v) be a solution of system (|])-(||) in the sense of Definition U^jl 



Thenv = Vf, where Vf is the function defined in (p5|), andu = dn inVtf = {x G 0,; Vf(x) > 
0}. 



The proof of Theorem 6A is essentially based on the techniques developed in [11, [|, 
O, ||. We will first prove the uniqueness of the first component of the solution of system 
(|l])-(2]). More precisely, we will show that if (u,v) is a solution of system ([l|)-(^), then 
u = dn in Qj := {x £ Q, : vj(x) > 0} (see Proposition ^4| below). 

Let us consider the functional $ : Hq (fJ) x L+(Q) — > R, where L 2 + {Q) = {z G L 2 (ft); z > 
0}, defined by 

(66) &(w,z) = — f(x)w(x) dx + z{x) [p(Dw(x)) — 1] dx. 

Jn Jn 

Lemma 6.2. If (u,v) is a solution of system @)-@, then (u,v) is a saddle point of<&, 
in the sense that 

z) < v) < $(w, v) V(w, z) G H%(Q) x L%(n). 
Proof. Since (u,v) is a solution of (||)-@, then 

v(x) [p(Du(x)) — 1] dx = 



and 

-2 



z(x) [p(Du{x)) - 1] dx < 0, \/zEL%{Q). 



Hence, for any z G Lq_(f2) we have 

v) = — f(x)u(x)dx 
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2 



(67) 

> - I f(x)u(x)dx + I z(x) [p(Du(x)) - 1] dx = <£(u, z). 



Moreover, by the convexity of p, for any w G Hq(Q) we have 

(68) p(Dw(x)) - p(Du(x)) > (Dp(Du(x)), Dw{x) - Du{x)) 

for a.e. i6fi. By Remark [4.2| , we can choose ip = w — u G Hq ($7) as test function in (| 
obtaining 

f(x)(w(x) — u(x)) dx + v(x) (Dp(Du(x)), Dw(x) — Du{x)) dx = 0. 



Thus, by ©, for any w G H^{Q), 

$>(w, v) — 3>(u, v) > — / f(x)(w(x) — u{x)) dx 
(69) J J? 

+ / u(a;) (Dp(Du(x)), Dw(x) - Du{x)) dx = 0. 



Collecting together ( p7|) and (p9|) we get the conclusion. □ 
In what follows we shall use the set of functions 

(70) LiPp(^) := {w G Lip(H); w = 0on dtt, p{Dw) < 1 a.e. in VL}. 

It can be checked that w G Lip^(Q) if and only if w G Lip(f2), w = 0on 30, and 

(71) w(x) — w(y) < p°(x — y) Vx, y G with [x, y] C f2 
(see |6|, Chap. 6]). 

Lemma 6.3. If (u,v) is a solution of system (|l])-(|2]), then the following hold. 

(i) u = dn in supp(/). 

(ii) (dn,v) is a solution of<^)~@). 

Proof, (i) By the maximality property of viscosity solutions we have that u < do, in Q. 
On the other hand, by Lemma |6.2| , &(u,v) < &(w,v) for any w G Lip*(f2). In addition, 
we have 



Q(u,v) = — f(x)u(x)dx, 
Jn 

Q(w,v) = — J f{x)w{x)dx+ I v(x) [p{Dw{x)) — 1] dx 
< — / f(x)w(x) dx . 



f(x)w(x) dx < / f(x)u(x)dx \fw G Lip^(O) . 



Then 



Choosing w = dn, we obtain that u = dn on supp(/). 
(ii) From (i) and Lemma |6.2| we have that 



$(dn,v) = $(u,v) < <f>(w,v). \/w G H^(n). 
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Hence, for every test function <p £ C^°($l) and every ft>0we have that 
< $(dn + hip,v)- <&(dn,v) 

= -h / f(x)tp(x) dx + / v(x)[p(D(dn(x) + htp(x)) - p(Ddn(x))] dx . 



< pplU \\D<p\\ LO o {n) 



Since p is convex we have 

p(D(d n (x) + h<p(x))) - p(Dd n {x)) 



h 

for a.e. x G and every h > 0, where 

pplU :=sup|£)p(OI = max |£>p(£)| 

due to the positive 0-homogeneity of Dp. Hence, from the differentiability of p in R n \ {0} 
and the dominated convergence theorem we get 



< - / f(x)ip(x)dx+ / v(x) (Dp(Ddn(x)), Dip{x)) dx . 
Jn Jn 

Replacing tp by —ip we also get the opposite inequality. □ 

Proposition 6.4. If (u,v) is a solution of system (|l|)-(0), then u = d^ in the set Qf = 
{i£!l; Vf(x) > 0}, where Vf is the function defined by (|35|). 

Proof. Let x € Sit Q f2 \ S. By the definition (35) of vj, and taking into account that 
M x (t) > for every t £ [0, r(x)) and / > 0, we deduce that there exists to £ (0, 7"(x)) such 
that, at the point xo = x + to Dp(Ddn(x)), one has /(xo) > 0. Since 

dn(x ) = dn(x) +t = da(x) + p°(x - x) , 

from this identity, Lemma |6.3| (i) , (fn]) and the inequality u < dn we get 

dn(x) = dn(x Q ) - p°(x - x) = u{xq) - p°(x - x) < u(x) < d n (x) , 
hence u(x) = dnix). □ 

Now that the uniqueness of the first component of the solution of system (||)-(|2|) is 
proved, it remains to prove the uniqueness of the second one. In order to do so, we will 
first exhibit for such a function a representation formula on the set Q \ £ and then analyze 
its behavior on S. 

Proposition 6.5. If (dn,v) is a solution of system (|])-(|2]), then for any zq G SI \ X and 

8 £ (0, t(zq)) we have 

v(zo) - v(z + 8 Dp(Dd n (z ))) M ZQ (8) 

(72) re 

= / f(z + tDp(Dd n (z )))M Zo (t)dt. 
Jo 

Proof. Let z\ = zq + 8 Dp(Ddn(zo)), let IL(zq) = {xo}, and define to = dn(zo) = p°(zo — 
xq), ti = dn(zi) = p°(z 1 - xo) = t + 8. Let Y : U -> W\ U C R n_1 open, a local 
parametrization of d£l in a neighborhood of xo, such that Y(0) = xq. Let $:WxK-> R n 
be the map 

(73) V(y,t)=Y(y)+tDp(v(Y(y))), (y,t) eU xR. 
Choose r > such that U(r) := {y £ IR n_1 ; \y\ < r} C U, and 

D(r) :={*(y,t); y € U(r), t £ Mi]} C fi \ £ . 
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The set D(r) can be viewed as a tubular neighborhood of the segment [2o>£i]- Let us 
define 

Si(r) = {*(y,ti); y G U(r)} , i = 0, 1, 
and let S^fV) denote the lateral surface of D(r), i.e. 

S 2 (r) = {V{y,t); y G dU(r), te[t ,h]}. 

All these surfaces are of class C 1 and are oriented with the outward normal with respect 
to D(r). 

For e > small enough let ij) e : R — > R and rj t : M n_1 — » R be the functions defined by 



if t < t or t > ti, 

1 if t G [t + e,ti - e], 



*=* iftG(t ,t + e ), if " y "^- e ' 



t_i=i if t e (ti — e,ti). 
Let be the function defined by 



if |y| > r, 
1 

7 ~\ y \ if r — e < \y\ < r. 



<p e (x) := 



ip e (t)r]e(y), if 3 y G I7(r) and t G [t , ti] s.t. x = (j/, t), 
0, otherwise. 



It is clear that <p e belongs to Lip(f2) and has support contained in D(r), hence can be used 
as test function in (|34]) . 

It is plain that <p e converges monotonically to 1 in the interior of D(r) as e — > + , hence 

(74) lim / f ip e dx = f dx . 

e^0+Jn JD(r) 

Let us compute the right-hand side of (|3^ ) when Lp = Lp e and u = dn- On -D(r) the test 
function <p e is defined by the relation 

<p £ (y(y,t))=><p e (t)ri e (y), y£U{r), te[t ,tx]. 

Differentiating the relation above with respect to t and recalling the definition ( |73|) of ^ 
we obtain 

(Dip e (*(y, t)), Dp(Dd n my, t)))) = ^(t) y G C/(r), t G [t , ti] . 

Then, taking into account that 4>' e (t) = for i G (to + e, ti — e), we get 

u (Dp(Ddn), DLp e ) dx 

n 

= [ [ 1 «(*) <Z)p(Ddn(*)), £Ve(*)> det DVdtdy 

(75) ^l/(r) -/to 

/ v(V)ip' e (t)r) e (y) detDVdtdy 

U(r) J t 

I (e) + h(e) + I 2 (e) + I 3 (e) , 
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where 

1 r rto+e 



J ( e ) = - f f° v (y) det D^dtdy, 
e JU(r-e) J t 

I x (e) = -- f f 1 v ($) det D^dtdy . 

e JUtr-e) JU-e 



£ JU(r-e) Jti-e 
1 f fto+£ 



h(e) = - f f° »(f)ii e (y) det dt dy , 

£ JU(r)\U(r-e) J to 

h(e) = --f f J v(V)r] e {y) detDVdtdy. 

e JU(r)\U(r-e) JU-e 



£ JU(r)\U(r-t) Jti-e 

Since t> and D^> are bounded in D(r), an explicit computation leads to \l2(e)+Is(e)\ < Ce. 
Passing to the limit in (ff^), by the continuity of v(^f) and we obtain 



lim f v (Dd Q , D<p e ) dx = VY-l)* / det D*(y, t;) dy . 

Recalling (frl), we finally obtain 

(76) / f(x)dx = J2(-iy f v^(y,U))detD^(y,U)dy. 

JD(r) f^O Ju ( r ) 

As a last step we want to pass to the limit as r — * + . From the continuity of v and 
we get 

(77) lim —-^ [ v(V(y, U)) det DV(y, U) dy = v{ Zi ) det (0, t;) , 

r^0+a{r) Ju(r) 

where cr(r) = u> n _ir n_1 is the area of the ball with radius r > in R n_1 . Finally 

lim — — f f(x)dx 
r^oo-(r) J D {r) 

lim — f f 1 /(*(?/, t)) det Dtf (y, t) 
r->0 air) Ju(r) Jtn 



r-*0 cj(r) 7l/(r) /to 
ti 

/(*(0,t)) det£>^(0,t)dt 

to 



(78) 



f(z + tDp{Ddn(z ))) det£>*(0,t + t)dt. 
From Lemma 4.10 in |10[] we have that 

det D*(0, t) = VGp(v(x )) det(J n _i - tF(i )), * G [to, *i] , 

where G is the determinant of the matrix of the metric coefficients. Collecting together 
@, ©, recalling the identity ©, and dividing by det 1^(0, t ) we obtain 

det(J n _i - hWixo)) 

v{zo) - v{zi) — — 

det(J n _i - t W{x Q )) 

f 6 f( ,. n m , , u , det(J n -i-(to + t)W(xo)) J4 

= / f(zo + tDp(Ddn{zo))) — — — — — (it. 

7o det(/ n _i - t W(xq)) 

The representation formula ( |72| ) now follows from ( p6| ) and the definition (32) of M ZQ . □ 
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For the proof of Proposition 6.7 below we need two more technical ingredients. The 
first one is the regularity result proved in Theorem |5.1| . The second one is the following 
convergence lemma due to H. Brezis (see ||, Theorem 1]). 

Lemma 6.6. Let 7: R n — ► R be a strictly convex function, satisfying the linear growth 
condition 

7(0 > co|e| - b , V£ G R n , 

for some positive constants bo and cq- Let (v,k)k C [-L 1 (f2)] n be a sequence of functions 
converging to u G [L 1 (f2)] n in the weak L 1 topology, and assume that j(u), j(uk) G 
[L 1 (r2)] n for every k G N. If lim^ Jq ^y(uk) dx = f^{u)dx, then (uk)k converges to u in 
the strong L 1 topology. 

Proposition 6.7. If(dn,v) is a solution of system (|l])-(|2|), thenv(x) = for every x G X. 

Proof. Since v is a continuous function, it suffices to prove that v = on S. Let us fix any 
xo G £ and choose e > sufficiently small such that B t (xo) C $7. Then, for any x G £?i(0) 
set 

d € (x) := dn ( x ° + ex ^ — da x ° t Ve ( x ) ■= v(x + ex), f e (x) := f(x + ex). 

By construction, for any e > as above d e (0) = and 

p(Dd e (x)) = p(Ddn(x + ex)) = 1 for a.e. x G B 1 (0) . 

Hence, there exist a sequence (e,),-, ej — > + and a Lipschitz function do : -£>i(0) — > R such 
that {d ej )j converges to do uniformly in -Bi(O). Moreover, since p(Dd ej (x)) = 1 in the 
viscosity sense in -Bi(O), by [Q, Proposition 2.2] also p(Ddo(x)) = 1 in the viscosity sense 
in -Bi(O), which gives p(Ddo(x)) = 1 almost everywhere. Since {Dd tj )j is bounded in L°°, 
we can also assume that it converges to Ddo in the weak L 1 topology. 
From in Lemma |5.3| , we have that the function 7^) = p(£) 2 /2 is strictly convex in 
R n . Since j(Dd ej (x)) = 1/2 and j(Dd (x)) = 1/2 for a.e. x G B^O), we have that 

lim / j(Dd e Jx)) dx = = / j(Ddo(x))dx. 

3^ooJ Bl (o) 2 7Bi(0) 



Recalling that {Dd ej )j converges to Ddo i n the weak L 1 topology, from Lemma 6.6 we 
conclude that (Dd t . )j converges to Ddo m the strong L 1 topology. 

Finally, the functions v t . and / e . defined above uniformly converge to v{xq) and /(xo) 
respectively and the pair (d €j , v tj ) solves 

(79) -div(v ej Dp(Dd ej )) = e j f ej in B 1 (0) 

in the sense of distributions, due to the fact that (dn,v) solves (§)-(§). Upon observing 
that Dp(Dd ej ) converges to Dp(Ddo) in the strong L 1 topology, we can pass to the limit 
as j — ► 00 in (|79|), obtaining that do is a weak solution of 

-div(v(xo)Dp(Ddo(x))) = x G #i(0). 

Now, if v(xo) 7^ 0, then do must be a solution to 

j-div(Dp(Dd )) = in Si(0), 
1 p(-Ddo) = 1 a.e. in Bi(0). 



From Theorem 5.1 we have that do G C 1,a (-Bi(0)). On the other hand, do cannot be 
differentiable in x = 0, because do is the 'blow up' of the distance function around a 
singular point xq. Hence v(xq) = and the proof is complete. □ 
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The last two propositions allow us to prove Theorem 6.1 as a simple corollary. Indeed, 
we already know by Proposition that if (u,v) is a solution of system (|l])-(^), then 
u = cIq on the set £lf = {vf > 0}. So it only remains to prove that v = v* in S7, where v* 



is given by (^). Proposition |6.7| guarantees that v = in S, while Proposition 3.5 implies 
that for any z$ E \ £ and E (0, t(zq)) 

v(z ) - v(z + 6 Dp(Dd n (z ))) M ZQ (9) 

f(z + t Dp(Dd Q (z ))) M Zo {t) dt . 



o 



Hence, letting 6 — > r(zo) and using the continuity of w we obtain that v(^o) = u/(zo)- 
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